Thermodynamics of ideal quantum gas with fractional statistics in T> dimensions 
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We present exact and explicit results for the thermodynamic properties (isochores, isotherms, 
isobars, response functions, velocity of sound) of a quantum gas in dimensions V > 1 and with 
fractional exclusion statistics < g < 1 connecting bosons (g = 0) and fermions (g — 1). In 
T> = 1 the results are equivalent to those of the Calogero- Sutherland model. Emphasis is given 
to the crossover between boson-like and fermion-like features, caused by aspects of the statistical 
interaction that mimic long-range attraction and short-range repulsion. The full isochoric heat 
capacity and the leading low-T term of the isobaric expansivity in D — 2 are independent of g. The 
onset of Bose-Einstein condensation along the isobar occurs at a nonzero transition temperature 
in all dimensions. The T-dependence of the velocity of sound is in simple relation to isochores 
and isobars. The effects of soft container walls are accounted for rigorously for the case of a pure 
power-law potential. 



I. INTRODUCTION 

A whole new arena for the study of effects of dimen- 
sionality has opened up in the wake of the advances 
in instrumentation achieved in the research that led to 
the observation of Bose-Einstein condensates [I], [2] . The 
magnetic and optical traps that were developed and per- 
fected in this and related lines of research can also be 
made in shapes that effectively constrain the kinematics 
of trapped gas atoms to fewer than T> = 3 dimensions 



In cigar-shaped traps, for example, the energy states 
of single particles with nonzero momentum components 
perpendicular to the axis are frozen out for the most part 
at low enough temperature. In disk-shaped traps the only 
energy states with significant occupancy at low T are 
those with zero momentum component perpendiclar to 
the plane. We are then dealing, effectively, with quantum 
gases in dimensions T> = 1 and T> = 2, respectively. A 
further effective change in spatial dimensionality can be 
roduced by controlling the firmness of the trap walls 



Of no less importance than effects of dimensionality 
are, of course, effects of interaction. They are notoriously 
difficult to handle in any effort that goes beyond mean- 
field theory or low-order perturbation calculations. Here 
a little explored non-perturbative approach presents itself 
as a promising alternative: the thermodynamic analysis 
of statistically interacting quantum gases. Statistical in- 
teraction grew out of the concept of fractional exclusion 
statistics, introduced as a tool to interpret the quasipar- 
ticle composition of the eigenstates in exactly solvable 
quantum many-body systems . 

The thermodynamics of statistically interacting de- 
grees of freedom is amenable to exact analysis under very 
general circumstances [13 . A most remarkable link ex- 
ists between the dynamical interaction in the form of a 
coupling term in the many-body Hamiltonian and the 
statistical interaction in the form of a generalized Pauli 
principle. For certain solvable models in T> — 1 one can 



replace the former by the latter and arrive at equivalent 
thermodynamic properties 

The chances for extending any known exact solution 
for a dynamically interacting quantum gas to T> > 1 are 
slim because the criteria for purely nondiffractive scat- 
tering |l2j are unlikely to be satisfied apart from highly 
contrived scenarios. By contrast, extending the exact 
solution for a statistically interacting quantum gas to 
higher dimensions is straightforward conceptually if not 
technically. 

Notwithstanding the absence of an exact correspon- 
dence between dynamical and statistical interactions in 
X> > 1, the exact thermodynamics of statistically inter- 
acting degrees of freedom has the potential of displaying 
a richness in phenomena unrivaled by mean-field theory 
or by perturbative approaches. Most importantly, it in- 
corporates a full and consistent account of fluctuations. 

Here we consider a statistical interaction limited to 
particles with identical momenta, in which case it re- 
duces to a statistical exclusion condition. The two pa- 
rameters of the system explored here are the coefficient 
of exclusion statistics, < g < 1, which spans a bridge 
between bosons (g — 0) and fermions (g = 1), and the di- 
mensionality T> > 1. While this system does exhibit the 
attributes typical for an ideal quantum gas, several of the 
highlighted features are indicative of specific aspects of 
the statistical interaction including features of long-range 
attraction and short-range repulsion. Selected results are 
known from previous work 0, ItH ITEl fl6| . 

A brief review of fractional exclusion statistics and its 
use in statistical mechanics is found in Sec.^ The appli- 
cation of these statistical mechanical tools to the model 
system under consideration is described in Sec. IIIII The 
exact thermodynamics of the ideal quantum gas with 
fractional exclusion statistics in a P-dimensional box 
with rigid walls is presented in Sec. IIVI The impact of 
soft container walls as might be relevant in the context 
of atomic traps is discussed in Sec. |Vj The main results 
are assessed in Sec. I VII as benchmarks for corresponding 
results pertaining to quantum gases with less constrained 
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statistical interactions as are forthcoming from work in 
progress [l7|. 

II. FRACTIONAL EXCLUSION STATISTICS 

The core relation of exclusion statistics is a generalized 
Pauli principle as introduced by Haldane [j| . It expresses 
how the number of states available to one particle is af- 
fected by the presence of other particles: 



(1) 



where the index i enumerates distinct particle species (a 
flexible notion) and the <?.; are coefficients of (fractional) 
exclusion statistics. For free bosons we have gt = and 
for free fermions gi — 1. Integrating Eq. Q yields the 
dimensionality of the one-particle Hilbert space in the 
presence of ni particles of species i: 



di = Ai - gi(rii - 1), 



(2) 



where Ai are constants. The dimensionality of the many- 
particle Hilbert space is 



W = 



n 



(3) 



where each factor represents the simple combinatorial 
problem of placing ni particles among di + ni — 1 dis- 
tinct states available to them. 

A system of (dynamically) free quasiparticlcs with frac- 
tional exclusion statistics has two principal specifications: 
(i) a set of orbitals at energies ti and (ii) a set of statisti- 
cal exclusion coefficients gi. The grand partition function 
is of the form [3 El 



Z 



n 



w, 



(4) 



where the Wi are determined by the nonlinear algebraic 
equations 



A 4 Wl . x f'l+Wi 
m(l + Wi) - gt In 



(5) 



k B T \ Wi 

The energy level occupancies (rij) are inferred from the 

Wi as 

, . 1 



(6) 
Wi + gt 

To illustrate the link between the combinatorial ex- 
pression and the statistical mechanical expression Q 
we consider an open system of M orbitals populated by 
(dynamically) free particles with (uniform) exclusion co- 
efficients g > 0. The number of distinct n-particle con- 
figurations implied by © is 

'A + (l-g)(n-iy 
n 



W n (M) = 



T(A + n-g(n-l)) 

r(n + r)r(A-g(n-i))' 



(7) 



where A is to be adjusted so as to yield integer values 
for W n (M) and a meaningful maximum capacity n max . 
For integer g we use A = M. For g = ^ we use A = 
M if M is odd and A = M + \ if M is even. This 
scheme is naturally extended to other fractional values of 
<7, limiting the maximum capacity to n max < (M + g)/g. 
The dimensionality of the Fock space, 



W tot (M)= ]T W n (M), 



(8) 



n=0 



can be expressed, for integer g, as a (higher-order) Lame 
sequence, 

W tot (M) = M+l (9) 

for M = 0, 1, . . . , g — 1 and 

Wtot (M) = Wtot {M — 1) + W tot (M - g) (10) 

for M = g, g + 1, . . . In the case of fermions (g = 1) we 
obtain W to t (M) = 2 M and for g = 2 a string of Fibonacci 
numbers: W to t(M) = Fm = 1,2,3,5,8,... For semions 
(.9 = \) we use 0- The result is a string of alternate 
Fibonacci numbers: Wt t{M) — F 2 m+i = 1,3,8,... 

The quantity we investigate for our illustration is the 
limit M — > oo of the proportion by which the number 
of many-particle states increases if we add one orbital to 
the system: R M = W to t{M + l)/W to t(M). For fermions 
we have Rm = 2 independent of M as expected. For 
semions we obtain the square of the golden section: 



Rm = 



F2M+3 M^oc 



1 



-(V5+1 



(11) 



F2M+1 

Now we derive R^ directly from Q and JSJ at T — 00: 



1 



W 



w 9 (l + w) 1 - 9 = 1. 



(12) 



The solution for fermions is w = 1, yielding R^ = 2 as 
before via combinatorics. In the semion case the asymp- 
totic result in 111|) is recovered from the solution of the 
quadratic equation in i|12|) . 



III. QUANTUM GAS WITH FRACTIONAL 
STATISTICS 

Before we go on to discuss how Eqs. Q © are applied 
to a quantum gas in T> dimensions, we turn to an exactly 
solved model in T> = 1. 



A. Calogero-Sutherland model 

The Calogero-Sutherland (CS) model describes mas- 
sive particles on a ring of circumference L with a (dy- 
namical) pair interaction (in units where h 2 /2m = 1) 
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o 2 



2g(g-i) 

d 2 



Here 



-sin (y(xi - Xj) 

7T \L 



(13) 



(14) 



represents the chord distance on the ring between par- 
ticles. Varying the parameter g across the interval < 
g < 1, amounts to a parametric link between free bosons 
(g = 0) and free fermions (g = 1). 

The spectrum and thermodynamics of the CS model 
can be described by an asymptotic coordinate Bethe 
ansatz |T^ . The thermodynamic Bethe ansatz inferred 
therefrom expresses the grand potential in the form 

fi(T, L, fi) = -k B T J +OC dk In (l + e-'W"° T ) , 

(15) 

where c(k) is determined by the Yang- Yang- type equa- 
tion m 



e(k) = k 2 - p 
k R T r+oc ' 



2tt 

with kernel 



/ dk'K(k~ k')\n(l 

J — oo 



e{k')/k B T 



K(k - k') = 2tt(1 - g)S(k - k'). 



(16) 



(17) 



The distribution of particles, pp(k), is found to be the 
solution, for given e(k), of the Lieb-Liniger-type equation 

m 



2irpp(k) 



1 + e e(k)/k B T 



= 1 



+ 00 



dk' K(k-k')p P (k'). 

(18) 



B. Generalization of CS model 

Returning to statistical interaction we consider a non- 
relativistic gas in a rigid box of dimensionality D and 
volume V = L v with energy- momentum relation eo(A;) = 
|k 2 | and exclusion statistics 

g(k-k')=gS(k-k'), (19) 

The grand potential derived from (0} becomes 

n = -k B T(f) D />fclni±^, (20) 
\2n J J w k 

where % is the solution of 

\k\ 2 -p 



k B T 



ln(l + «fc)-0lni±^. (21) 



The particle density in k-space is represented by the func- 
tion 



("k> 



w k + g 



(22) 



The fundamental thermodynamic relations, from which 
most thermodynamic properties are conveniently de- 
rived, are integrals involving Wk and (rik): 



k R T V27T 



v 



PV f L ^" fpkJnl + WL, (23) 



TV : 



P 



dPk (rik) 



1 l^)"/' k • 



(24) 



(25) 



These relations state the dependence of the pressure p, 
the average number of particles Af, and the internal en- 
ergy U on fugacity z, temperature T, and volume V. 
They are also known as the thermodynamic equation of 
state (H2J)-(|2U) and the caloric equation of state (|2"4*|) - (|2*5|l . 

It was recognized 0, 0, that the thermodynamic 
Bethe ansatz solution l|15 |) -l(T5 )l of the CS model is equiv- 
alent to the solution of a statistical interacting gas as 
presented above for T> = 1 if we make the following iden- 
tifications: 



w k = exp 



\k B T 



(n k ) = 2irp P (k), 



(26) 



g(k,k') =5(k-k') K(k-k'). (27) 

2n 

For the thermodynamic extension to T> > 1 of the CS 
model we use the statistical interaction (|19|) and the den- 
sity of 1-particle states of the nonrelativistic ideal gas, 



L\ V TT^ 2 v/2 _ 1 



D ° {£a) ~{2^J T{V/2j^ 
to rewrite Eqs. (20)1 - (|2*2")l in the form 

\l = -k B T deo£>o(eo)ln , , , 

Jo w(e ) 

iexp(^) =KeoF[l + ^ )] 1 - 9 , 
z = cxp(/z/fc s T), 



w(eo) + .9' 



(28) 

(29) 

(30a) 
(30b) 

(31) 



and to reduce Eqs. I|23|l - H25(l to integrals over functions 
of e . 
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C. Average level occupancy 



In ideal quantum gases the average level occupancy 
(n(eo)) is a unique function of (eo — n)/kBT independent 
of T>. A plot of this function for the generalized CS model 
is shown in Fig. ^ for various values of g. In the boson 
case we have iu(eo) = e^ e °~ p ' ) / fcBT — 1 and in the fcrmion 
case w(eo) = e^°~^/ kBT . For semions we have 



V 



= G v/2 (z,g) 



V 



V 

2 kBT 



Gv/2+i(z,g), 



(36) 



(37) 



w(e ) 



\ e (eo-^)/k B T] 2 . i _ i 
L J 4 2 



(32) 



where 



The exact analytic solution of (|30H with < g < 1 is 
known in a series representation [16j. The asymptotic 
level occupancy for (eo — ^/ksT — > — oo is l/g. 




(£ -(I)/k B T 

FIG. 1: Average level occupancy (n(eo)) versus (eo -/iJ/feT 
for 3 = 0, 0.1, 0.25, 0.5, 0.75, 1 (from top down). 



IV. THERMODYNAMICS OF THE 
GENERALIZED CS MODEL 

Exact and explicit results for the thermodynamics of 
the generalized CS model can now be calculated from the 
general expressions derived in Sec. 11111 



A. CS functions 



Introducing the CS functions 



G n (z,g) = 



dx 



r(n) J w(x)+g' 



[wWll + wix)] 1 -" = — , 

z 

for n > we rewrite Eqs. I|23|) (1251) as 



(33) 
(34) 



= G v/ 2+i{z,g), 



h 2 



S 2 /2m=l 



27rrnfcsT 



47T 



(38) 



is the thermal wavelength and the term in (|36() enclosed 
by brackets is relevant only if g = and V > 2. For 
9 = and g — 1 the CS functions become the Bose- 
Einstein (BE) functions and Fermi-Dirac (FD) functions, 
respectively. The range of fugacity is < z < 1 for 
bosons. For all other cases z has no upper bound. 

The overall shape of the CS functions is illustrated in 
Fig. |21 For \ < g < 1 the curves are concave. For 
< g < g they have a convex portion at small z and 
then switch to to concave behavior. The curve for g = 
is convex over its entire (restricted) domain. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



FIG. 2: CS functions @3 for 



(35) g = 0, 0.1, 0.25, 0.5, 0.75 (from top down). 



plotted versus z for 
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The CS functions have the following power series ex- 
pansions (lt| : 

r , x f r(Z - iff) 



The familiar recurrence relation for FD and BE functions 
is valid for the CS functions in general: 



£j i» r(Z)r(i - iff) 



9 

ZT^G„ + i(z,sr) = G n (z,g). 



(42) 



= z + |^(l-2ff) + |^(l-3ff) ) CO, : ! ,. 



The radius of convergence depends on g: 

1 



-(5) = 



9 S(l- fl )l-9- 



(39) 



(40) 



The first two terms of the asymptotic expansion for large 
z at g > are jlHll^ 



G n {z,g) 



(lnz) n 
gT(n+l) 



it 2 gn(n — 1) 
6 (lnz) 2 



(41) 



We also use l|42[) to extend the definition l|33[l to n < 0. 



B. Equation of state 

The functional relationship between pV/NksT and fu- 
jacity z, 



pV G v/2+ i(z,g) 



Nk B T G v/2 (z,g) ' 



(43) 





FIG. 3: Equation of state 03) in V = 1, 2, 3, 5 for g = 0, 0.01, 
0.1, 0.25, 0.5, 0.75, 1 (from bottom up). Note the different 
vertical scales. 



as inferred from l|35|) and H36(l . highlights the deviations suits of pV/JVksT versus z are presented in Fig. 03 In- 
from Maxwell-Boltzmann (MB) behavior. Numerical re- 



6 



creasing z means, for example, decreasing T at fixed 
Af/V or increasing Af/V at fixed T. Any downward (up- 
ward) deviation from pV/AfksT = 1 is suggestive of a 
boson-like (fermion-like) feature. 

The curves for g > \ exhibit a fermion-like monotonic 
increase over the entire range of z. The curves for < 
g < i start out with (boson-like) negative slope. They 
reach a smooth minimum and then grow without bound 
in fermion-like manner. The case g — is exceptional. 
Here the curve ends in a cusp singularity at z = 1. In 
V = 1 and V = 2 the pressure is zero at the singularity. 
The endpoint of the curve occurs at T = or Af/V = oo. 
In T> > 2 the pressure is nonzero at the singularity, and 
the endpoint of the curve signals the onset of a Bose- 
Einstein condensate (BEC). For fixed Af/V this happens 
at T > 0. The condensation proceeds along the dashed 
line. 

The physics underlying the crossover between boson- 
like and fermion-like behavior may be interpreted by at- 
tributing to the statistical interaction of the generalized 
CS model a long-range attractive part and a shorter- 
range repulsive part. The repulsive core is present for 
all g > 0. The attractive part is only perceptible for 
g < i. The former prevents the collapse of the system 
into a BEC at g > 0. The latter causes a negative in- 
teraction pressure for g < | (a reduction of p relative 
to the kinetic pressure of the MB gas) if the average in- 
terparticle distance in units of the thermal wavelength is 
sufficiently large. 

The interplay between the two parts of the statisti- 
cal interaction upon variation of the parameters g and 
T> produces a host of interesting thermodynamic effects. 
Their appearance in isochores, isotherms, isobars, re- 
sponse functions, and in the velocity of sound will be 
discussed in the following. For this purpose we intro- 
duce reference values for temperature T, reduced volume 
v = V/Af, and pressure p, based on the equation of state, 
pv = ksT of the MB gas and the criterion that the av- 
erage reference volume per particle is a hypercubic box 
with sides equal to the thermal wavelength (|38l) . For 
isochoric, isothermal, and isobaric processes we thus use 



Isochores 



The dependence of p on T at constant v can be ex- 
tracted from (|35|) and i|36|) : 



P _ Gv/2+i{z,g) 



Pv [G v/2 (z,g)Y+W' 
T 1 



T v lG v/2 (z,g)]yv 



(44a) 
(44b) 



with no restriction on z for < g < 1. In Fig. 01 we 
show isochores for various g and D. Each curve repre- 
sents a universal isochore, valid for arbitrary values of v. 
The shape of the curves at high T reflects the emerging 
MB behavior. The leading correction term in the high-T 
expansion of (|44H . 



P T- 
Pv 



T 



1/2 -g (Tv 

2 V/2 [ rp 



V/2 



(45) 



describes bundles of isochores whose vertical separations 
increase with increasing T/T v if T> < 2, stay constant 
if I? = 2, and decrease if I? > 2. The deviation from 
the MB isochore, p/p v = T/T v , is negative for g < | 
and positive for g > decreasing in magnitude as T> 
increases. For semions the deviation is of higher order in 
T v /T. 

The intercept of the isochore at T = as extracted 
from the leading term in the asymptotic expansion of 
(H2J is 



Po 

Pv 
f, 



T v V/2+1 



^ = \gY{V/2 + l)] 



2/V 



(46a) 
(46b) 



where /iq — ksT v is the chemical potential at T = 0. If 
we replot the data for g > in Fig. 0] as p/p v versus 
T/T v with p v = p v (T/T v ) we obtain the MB isochore in 
the limit T> — > oo. 

In the boson case, Eq. (|44|l is still valid for z < 1. This 
covers the full T-range in T> < 2. In T> > 2 in the limit 
z —* 1 occurs at 



Vt 



k B T x 
-in 



4?r 



V 2/V ' Pv V 2/V+1 
T>/2 , kBT)V /2+l 
, PT ~ 



4tt 



k B T ) ' " J (4ir) v / 2 
ksTp = (Att)^p^ , v p = (4tt/p)^ 



(v = const.) 

(T = const.) 
(p = const.) 



These reference values are well-behaved in the boson 
limit, which is an advantage for our comparative plots. 
For some purposes (limited to g > 0) an alternative 
choice, based on the chemical potential at T = 0, will 
be more convenient. 



Pc CP/2 + 1) 
Pv [C(2?/2)P/f+i' 

Tc = 1 
T v [C(P/2)] 2 / 1 " 



(47a) 
(47b) 



where G n (l, 0) = £(n) is the Riemann zeta function. The 
values p c , T c mark the onset of a BEC. Note that T c be- 
gins to rise from zero at T> = 2 and approaches the value 
T v as T> — > oo. The bosonic isochore in the coexistence 
region is 



P_ 

Pv 



V/2+1 



((V/2 + 1) (T < T c ). 



(48) 
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The singularity at T c > becomes stronger as V in- 
creases. For 2 < T> < 4, there is a discontinuity in curva- 



ture. It turns into a discontinuity in slope for V > 4 and 
becomes a discontinuity in the function itself for T> = oo. 



J 




°- 0.6 




FIG. 4: Isochores in D = 1,2,3,5 for g = 0, 0.01, 0.1, 0.25, 
0.5, 0.75, 1 (from bottom up). In T> = 1,2 the curves for 
g — 0, 0.01 are unresolved. Note the different vertical scales. 



r 



D. Isotherms 



The dependence of p on v at constant T in a parametric 
representation follows again from l|35|) and (|36|l : 



PT 



V 

Vt 



Gv/2(z,g) 



(49a) 
(49b) 



applicable for < g < 1 with unrestricted z and for g = 
with z < 1. In Fig. we show isotherms for various g 
and T>. The convergence of the curves in the low-density 
regime reflects the emerging MB behavior. The leading 
correction to Boyle's law, 



p v ^oo vt 
pr v 



_ 1/2 - g /Vr\ 
2 C / 2 V V J 



(50) 



is negative for g < \ and positive for g > i. It weakens 
in magnitude as T> increases. The divergence of the iso- 
chores for g > in the high-density limit is of the form 



P_ 

PT 



v-r 



(T>+2)/T> 



2/P [T(V/2 + I)]*/* 
5 2?/2 + l ' 15 j 



consistent with MB behavior for g > in the limit T> — > 
oo, if we switch to alternative reference values based on 
the chemical potential. 

For bosons the pressure p/pt remains finite as v/vt 
0. In V < 2 the limit z — ► 1 implies 



(52) 
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whereas in V > 2 at z = 1 we have 



v T CP/2) ' PT 



CP/2 + 1). 



(53) 



The intercept of the bosonic isotherm decreases with in- 
creasing T> and approaches unity for T> — > oo. In 2? < 2 



the bosonic isotherm is a smooth curve. Its slope at 
v/vt — is negative in T> < 2 and zero in 2? = 2. In 
V > 2 the pressure is constant at p c for < v < v c 
along the isotherm. In the limit T> — > oo we have the 
MB isotherm, p/pr — vt/v at v > vt joined by a the 
constant p/pt = 1 at v/vx < !■ 





E. Isobars 



T>. The leading correction to MB behavior at high T, 



The dependence of v on T at constant p is to be cal- 
culated from 



v [G v/2+1 {z,g)r/tv + 2) 



1 V 

T 
if 



Gv/2(z,g) 
1 



p [Gv/2+i{z,g)] 



2/(X>+2) ' 



(54a) 
(54b) 



applicable for < g < 1 with unrestricted z and for g = 
with z < 1. In Fig. we show isobars for various g and 



V_ T-joo _T 



1/2 -g /%\ 
2 p/2 y t / 



25/2+1' 



(55) 



explains the observation that the vertical separations of 
the curves in Fig. [5] decrease with increasing T/T p in all 
T>. Again the deviation from MB behavior switches sign 
at g = h . The intercept at T = of the isobars for g > 
is 



lim 

t^o v, 



fJ 



2/CD+2 



p 



[r(p/2 + i)] 2 /( p+2 > 

(Z>/2 + if7p5+3 ' 



(56) 
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again consistent with MB behavior for g > in the limit 
T> — > oo provided we use alternative reference values 
based on the chemical potential. 

The isobaric curves for bosons end in a critical point 

at 



Vc 

v p 

Tr 



{((V/2 + l)f/^ 



CP/2) 
1 



(57a) 
(57b) 



T P [CCD/2 + l)] 2 /^+ 2 ) ' 
Note that the critical (reduced) volume is nonzero only 



in T> > 2 whereas the critical temperature is nonzero 
in all T> > 1. At T < T c the boson gas is unable to 
maintain the prescribed pressure. In D < 2 the bosonic 
isobar terminates in a cusp at v/v p = 0. The switch 
from a pure gas at T > T c to a pure BEC at T < T c 
occurs when v/v p = 0. There is no coexistence region. 
In T> > 2 the bosonic isobar terminates in a cusp at 
T = T c > and v = v c > 0. Here the two phases 
do coexist. During condensation v/v p gradually goes to 
zero along the vertical dashed line. With increasing D 
both Vc/vp and T c /T p increase toward unity. 



F. Isochoric heat capacity 



Af 1 (dU/dT) v , derived from that expression reads 



From (|36|l and (|37(l we know that the scaled internal 
energy U /AfkgT v differs from the isochore (|44|l by a mere 
overall factor T>/2. The isochoric heat capacity, C v = 



k B 



V 
~2 



Gy/2+\{z,g) _ V G v/2 {z,g) 
G v/2 (z,g) 2 G. 



27/2-1 



(58) 

In Fig. [7| we show plots of C v /kB versus T/T v from 
(I44|l for various g and V. In T> = 2 the heat capacity 
is independent of g, a result known since 1964 for the 
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O 0.3 - 




FIG. 7: Heat capacity in V = 1, 2, . . . , 5 for g = 0, 0.01, 0.1, 
0.25, 0.5, 0.75, 1 (from top down or bottom up). In T> = 1 
the curves for g — 0, 0.01 are unresolved. In T> = 2 the heat 
capacity is independent of 3. In 7? = 5 the heat capacity 
has a discontinuity AC„/&b ~ 3.21 at T c . Note the different 
vertical scales. 



FD and BE gases and recently extended to < g < 1 



C v 7r 2 T T v y, / , 2 T , 2 T 2 



fc s 3 T„ T 



n T v n 2 T 2 



-nT„/T 








9 = / 




/ 0.01 
















D = 5 



0.2 



0.4 



0.6 
T/T„ 



(59) 



r 



In D < 2 only the curves for <? > i have a local max- 
imum and in T> > 2 only the curves for g < i. This 
double switch is reflected in the leading correction to the 
MB behavior at high T: 



C a 



T^oo 



V 

2" 



(1/2 - 5 ) (25/2-1) / T, 



2 V/2 



T 



V/2 



(60) 



The leading low-T behavior of C v is linear for all g > 0: 
C K t^o tt 2 2? 5 (^- 2 )/ c T 



6 [r(2?/2 + l)] 2 / p T t ,' 



(61) 



In T> < 2 (I? > 2), the initial slope increases (decreases) 
with decreasing g, in agreement with the numerical re- 
sults of Fig. □ 
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The well-known power-law behavior of the bosonic 
heat capacity at low T reads 



C v t^o V (V 
k B 2 I 2 ^ 



1 



27/2 



(62) 



V > 2 shows a qualitative change of behavior in V = 4. 
In 2 < 2? < 4 the heat capacity CV, has a local maximum 
at T c with a discontinuous slope. In T> > 4 the function 
itself becomes discontinuous at T c as is evident in Fig.0 
The size of the discontinuity is 



It represents the leading singularity of the exact result in 
T> < 2 and, at the same time, the exact result itself for 
< T < T c in T> > 2. The latter result can be rewritten 
in the form 



k B 2 V 2 



C(g/2 + l) 
C(2?/2) 



27/2 



(63) 



for < T < T c with T c from gTJ. 

Inspection of the bosonic heat capacity results l|58l) and 
(|63|l for T —t T c from above and below, respectively, in 



AC V _ V 2 C (Z>/2) 



fc s 4 C (2?/2 - 1) 



(V > 4). 



(64) 



There is no latent heat. Only in the limit T> — + oo a first- 
order transition emerges. In large T> the heat capacity 
stays very close to zero from T = up to the near vicinity 
of T c , where it shoots up to a high value of 0(T> 2 ). At 
T c it drops back to a value of 0(1?) and remains nearly 
constant. In the limit T> — + oo this spike turns into a 
delta function and its weight determines the latent heat. 




G. Isothermal compressibility with v/vt from (|49|) . In Fig.[S]we plot ptKt versus v/vt 

for various g and V. 

We obtain the following parametric representation of 
k t = —v~ 1 (dv/dp)T from (|49(l : 

v G v/2 -i(z,g) 

PTKT = J\ -, s-, 65) 
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The leading correction to MB behavior at low density, 
v/v T > 1, 



v — >oo V 

vt 



1 



1/2 - g V T 



2©/2-l 



(66) 



is strongest in low D. The compressibility curves for 
g > | display fcrmion-like behavior across the entire 
range of v/vt- The value of kt stays below the MB 
value and decreases monotonically with decreasing v/vt, 
approaching zero in the limit v — > 0, reflecting the re- 
pulsive core of the statistical interaction. For g < i by 
contrast, the curves start out above the MB line. Here 
the long-range attractive part of the statistical interac- 
tion is dominant, producing boson-like behavior. As the 
particle density increases a crossover from boson-like be- 
havior at low density to fermion-likc behavior at high 
density manifests itself as a shoulder or as a precipitous 
drop after a smooth local maximum. The leading high- 
density term of all compressibility curves for g > is a 



power-law with Independent exponent: 



V 



2gV' D [T(V/2 + l)} 2 /' D \v T 



(67) 



The uppermost curve in each panel of Fig. [8] is for 
bosons. It shares with all the other curves a decreasing 
initial trend as v/vt is lowered in the low-density regime. 
Then it goes through a smooth minimum (a property 
shared with curves for g « 1) and rises to either a di- 
vergence or a cusp at v c . In T> < 2 we have v c = 0. 
Here the isothermal compressibility exhibits a power-law 
divergence, kt ~ (v/vt)^ v ^ 2 ^ t '^ , in 1 < T> < 2, and an 
exponential divergence, kt ~ {v / vt) 2 e v ^ VT , in T> = 2. In 
2 < T> < 4 we have v c > as given in (|53|l and still 
diverges. However, in T> > 4 the divergence turns into a 



cusp at Ptc k Tc 
have kt = oo. 



C(2?/2 - I)/ l((V/2)] z . Atv<v c 



we 
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H. Isobaric expansivity 



We calculate a p = v 1 (dv/dT) p from the general re- 
lation a p — KT(dp/dT) v and the CS-specific relation 
C v = (V/2)v(dp/dT) v to arrive at 



rp _ ^J>_ 

lpCtp — 



G 



15/2+1 



(z,g)G- 



G-D/i{z,g)G v/2 {zig) 



p i -i{z,g) _ v 

2 

(68) 



with T/Tp from l|54|l . In Fig. [5] we show expansivity 
curves T p a p versus T/T p for various g and V. We ob- 
serve that the correction to MB behavior at high T is 
such that the expansivity is suppressed for g > i and 
enhanced for g < A: 



TpOtp 



T, 



(l/2- fl )(2?/2 + l) (T % 



2 -D/2 



V/2+l 



(69) 

The characteristic features of all compressibility curves 
for < g < 1 are a smooth maximum and a linear ap- 
proach to zero in the low-T limit. The maximum is flat 
near the fermion limit and becomes increasingly high and 
narrow as the boson limit is approached. For g <C 1 the 
position of the maximum is close to the bosonic T c as 
given in (|57|) . The slope of the linear low-T behavior of 
a p depends on g and V: 



Tpdp 



(D-2)/(U+2) 



6 [r(X>/2 + l)]4/(x>+2) Tp 



(70) 



Note the absence of any ^-dependence in T> — 2 and the 
opposite trends regarding ^-dependence in T> < 2 and 
T> > 2 reminiscent of trends seen in the heat capacity. 

The bosonic expansivity increases monotonically with 
decreasing T and ends in a singularity as T c is approached 
from above. In T> < 4 the singularity is a divergence and 
in V > 4 it is a cusp. 



I. Velocity of sound 

We start from the relation c = (/?ks) -1 / 2 for the veloc- 
ity of sound, where p — m/v is the mass density and n$ 
the adiabatic compressibility. We use standard relations 
between response functions to arrive at the following ex- 
pression, all in terms of dimensionless entities previously 
determined 1301: 



mc 2 (v/vt) 
k B T (j>tKt) 



{T/Tp)\v/v T ){Tpapf 

{PTK T ){C v /k B ) 



(71) 



For the CS model the right-hand side of l|71|) can be 
greatly simplified: 



G 



T>/2+l 



{z,g) 



Gv/2(z,g) 



7=1 + -. (72) 



How does the velocity of sound c vary with temperature 
T if we keep the (average) particle density 1/v or the 
(average) pressure p fixed? Remarkably, we find very 
simple universal relations between c and the isochore or 
isobar itself for the two situations, respectively: 



mc p 

——=7— (v = const.), 
k B ± v Pv 



—— = 7— (p = const.), 

KBlp V p 



(73) 



(74) 



with the dependence of p/p v on T/T v as discussed in 
Sec. IIV CJ1 and the dependence of v/v p on T/T p as dis- 
cussed in Sec. IIV El 

We conclude that c has a monotonic dependence on T 
no matter whether we keep v or p constant, for g > 
the velocity of sound stays nonzero in the limit T — > 0. 
In the boson gas c is affected by the onset of BEC. There 
are no sound waves in the condensate. The isochores at 
T < T c in Fig. 0| or the vertical portions of the isobars at 
T = T c in Fig. [S]give us information about c in the gas 
coexisting with the condensate. 



V. EFFECTS OF SOFT TRAP WALLS 

Consider an ideal quantum gas in T> > 1, trapped by 
a spatially isotropic power-law potential 0, \ 



I r \ r i 



(75) 



where r is the radial coordinate in 2?-dimensional space. 
The size of the trap is determined by the "width" R and 
the "depth" Uq of the potential well. The softness of the 
confinement is controlled by the exponent g. Harmonic 
traps have g = 2. Lowering the value of g makes the sur- 
rounding wall softer, increasing it makes the wall harder. 
For g — > oo, the wall becomes rigid and the bottom of 
the potential becomes flat, which corresponds to the sit- 
uation considered in Sec. IIVI 

Softening the trap walls affects the density of energy 
levels. For the power-law potential (|75|) expression 128|) 
must be replaced by 



D(e ) 
with 
V 



V 



2 v tt v / 2 T{V/2) 



D/2-1 
£ 



Qv(g) (76) 



r P/2 



r(x>/2 + 1) 



R L 



Qv{g) 



T{V/g + l)r(Z>/2) 
T{V/g + V/2) 



(77) 

In soft-wall traps it makes sense to consider processes at 
V = const instead of isochoric processes. Under these 
circumstances the gas expands when heated up. In gen- 
eralization of relations (|36|l and 1)3 7|) we now have 



Af=^T(V/g 



■D/r, 



G 



V/2+T>/ 



n(z,9), (78) 
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U 



Vk B T ( k B T 



x r 



v/r, 
V 



+ ~2 j Gc /2+^/'?+i( 2 '3)- ( 79 ) 



As before, an additive term 2/(1 — z) has to be considered 
in l|78l) if g = 0. We now use the reference temperature 



-in 



-P/2 



■£Jn 



x [r(0/r? + l)]W^+W5 



(80) 



For bosons the onset of BEC occurs at T c > in 
2? > 2? m = 277/(2 + 77). The marginal dimensionality 
decreases from 2? m = 2 to D m = 1 as the rigid-wall 
container softens into a harmonic trap. The transition 
temperature in units of T v is 



^ = [((D/t,+P/2)]'™ . (81) 



VI. CONCLUSION AND OUTLOOK 

We have explored the thermodynamics of the general- 
ized CS model in T> > 1 dimensions. In this model the 
statistical interaction is limited to pairs of particles with 
identical momenta. This reduces the coupling, effectively, 
to a statistical exclusion condition. Several observed phe- 
nomena suggest the presence of a long-range attractive 
part and a short-range repulsive part in this particular 
statistical interaction. 

The generalized CS model is found to preserve several 
attributes that are characteristic of the familiar ideal BE 



and FD gases - attributes which may very well count 
as hallmarks of ideal quantum gases in general: (i) The 
average level occupancy (n(eo)) is a unique function of 
(eo — fi)/kgT for given g, independent of T>. (ii) The fun- 
damental thermodynamic relations l|35l) - (|37[) for given g 
and with the prescribed structure on the left-hand side 
are unique functions of the fugacity z. (iii) The right- 
hand sides of Eqs. (|3*3|l and l|3"7|) are identical. 

The consequences include that the isochore and the 
(properly scaled) internal energy have the same depen- 
dence on T, that the quantity pV/J\fk B T has a unique 
dependence on z for given V and that the square of the 
(properly scaled) velocity of sound at constant v (con- 
stant p) has the same T-dependence as the isochore (iso- 
bar). The thermodynamics of the generalized CS model 
as described in Sec. lIVl is thus ideally suited in the role of 
benchmark for the thermodynamics of gases with more 
generic statistical or dynamical interactions. From work 
in progress on several model systems in T> dimensions 
with statistical interactions not restricted as in |1 we 
have compelling evidence that none of the ideal quantum 
gas hallmarks (i)-(iii) are upheld anymore. The stage 
is thus set for the exact analysis of qualitatively new 
thermodynamic phenomena including phase transitions 
in non-bosonic quantum gases. 
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